Mode interaction in multi-mode optical fibers with Kerr effect 
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We generalize the projection to orthogonal function basis (including polarization modes) method 
for nonlinear (Kerr medium) fiber and use this method in a case of two-mode waveguide. We 
consider orthogonal Bessel functions basis that fit the choice of cylindrical geometry of a fiber. The 
coupled nonlinear Schrodinger equations (CNLS) are derived. Analytic expressions and numerical 
results for coupling coefficients are given; the fiber parameters dependence is illustrated by plots. 

PACS numbers: 42.65.-k, 42.65. Tg, 42.81.-i, 42.81.Qb 



I. INTRODUCTION 

It is well known that Nonlinear Shrodinger (NS) equa- 
tion, similar to celebrated KdV, presents the most uni- 
versal physical embedding of soliton behavior. How- 
ever, its experimental realizations are still rather poor 
PJ and rather demonstrate the solitons than introduce 
a true physical tool for nonlinear phenomena investiga- 
tions. The main reason for this is its uni-dimensionality 
and uni-directionality - the embedding needs some kind 
of projecting with exact formal and experimental speci- 
fication. 

The important example of the theory is the optical 
guide (fiber) propagation that isolates modes due to 
boundary conditions in a fiber cross-section and the di- 
rection of propagation is fixed by a way in which the 
waveguide is excited. So the model based on NS is es- 
sentially one dimensional, hence the fiber propagation is 
perhaps, the best one for physical realization of soliton 
behavior of electromagnetic wave-train. Nevertheless the 
existing theory adopts the model of waveguide prop- 
agation, that looks rather "technical" than physical one. 
The authors use the parameter of NS equation as em- 
pirical ones without explicit expressions dependent on a 
whole geometry and physics of a waveguide. 

In the papers 0] and authors investigate polariza- 
tion behavior of light in a Kerr fiber; additionally the 
authors of Q took into account power exchange between 
modes. The authors of these articles declare the complete 
basis expansion for the problem solution but in fact use 
simplifications choosing for the basis functions (for higher 
order modes) the Gaussian functions. This approxima- 
tion is used for nonlinear coefficients (in our notation ©i) 
evaluation. The formulas for the coupling constants (cjj 
in p|) exhibit the fiber radius dependence (proportional 
to 1/Vq) that essentially differs from ours. 

By the work we started the detailed investigations; 
we described and discussed a method to derive CNLS 
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(Coupled Nonlinear Schrodinger) equations for a multi- 
mode fibers focusing on a mono-mode case (which leads 
to NS equations) and made numerical calculations for it. 
We elaborated an approach which do not apply poten- 
tials, each component of electromagnetic field in cylindri- 
cal coordinates is represented directly as a series in corre- 
spondent Bessel functions and exponentials Ji(ai n r)e lllfi . 
Such scheme introduces natural mode notion in a fiber 
cross-section which have obvious link with the standard 
one (TE,TM modes). 

Here we derive CNLS equations for two modes (where 
I = 0, ±1 and n — in both cases). Information about 
the parameters allows to derive the explicit dependence 
of soliton properties on geometry and material constant 
of a medium. We use the model of step index waveguide 
with different refraction indexes (without assumption of 
weakly guided fiber like in 0, 0, 0)- The main result 
of this part of our study is the analytical expression for 
mode interaction coefficients that exhibit rather strong 
dependencies on waveguide radius. 

We take the fundamental mode (which is known from 
linear a theory as HE11) considering it and the second 
mode (TE01) as vectors of our basis. Both modes have 
different group velocities, this difference influence on the 
coupling coefficients magnitude. This model could be 
used in other waveguides with different geometries (e.g. 
the elliptical one |5j)- Both single- mode fibers with two 
polarizations and multi-mode fibers are used in optical- 
switching devices ||. 

Considering two modes with different propagation con- 
stants (this mean 5^0), we expected that interaction 
between modes occurs only at the beginning of propa- 
gation (we inquire the propagation to begin from exci- 
tation at one of waveguide's ends). This can appear in 
the strong birefringent photonic crystal fibers (PCF). In 
a consequence this interaction can lead to change power 
intensity of modes. For example, such properties are very 
important in case of photonic crystal fibers (PCF) or dis- 
persion managed solitons (DM). 

It is necessary to point out that in the cylindrical 
waveguide even three modes can interact: HE11, TE01 
and TM01. Modes TE01 and TM01 have the same cut- 
off frequency and propagation constant . This case is 
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very similar to single mode with two polarizations case 
and we considered this in 0. In case of HE 11 - TE01 
mode interaction, it is needed to take into account two 
mode interaction with different propagation constants. 
In this work we investigate only this (HE11-TE01) inter- 
actions to simplify equations. 

We start in the second section from basic equations 
of linear and some of nonlinear theory of electromagnetic 
fields. In third section analytical solution of a given prob- 
lem is presented. In fourth section some of numerical 
calculations for nonlinear coefficients as a function of the 
guide radius are given. The last section is a compari- 
son with other cases of analyzed equations and general 
conclusion is presented. 



II. BASIC EQUATIONS 

Consider the propagation of optical pulse at isotropic 
medium and choose a cylindrical dielectric waveguide 
with small Kerr nonlinearity. The electric field can be 
written as 



Ei = i^e*" + c.c, 



(1) 



where i = x,y,z. We Introduce a linearly polarized field 

as HE3 



and the polarization vector component as 



(2) 



1-r 



\A m \* + ±(\A m \* + \A v \*) 



A z e iwt 



-A z (A 2 . + A 2 



(3) 



where A z is the complex conjugate of A z . 

The electric field component in Bessel function basis, 
which is standard for cylindrical waveguide can be writ- 
ten as pHH^ I 



2 ^— ' ai n 



(4a) 



Ey(x,y,z,t) 



2 7^ ain 

(4b) 

here ai n is eigenvalue for linear cylindrical waveguide 
problem 

a 2 = uj 2 e ^ ei - k 2 , r<r , (5a) 
(3 2 = k 2 - ui 2 e fi e2, r > r , (5b) 



where rg is waveguide radius and the variable amplitude 
A depends on propagation axis and time. 

Using the multi-mode model 0| (where the modes or- 
thogonality over the fiber cross section is used) we obtain 



(a z + a 2 Q1 )A p 01 



2eo/io 

7TiV 01 



ro 27r 



rJ Q (a Q1 r)—P z dipdr, 



(6) 



(D z + a 2 n )A p 11 



ttJVu 



r 2tt 



rJ 1 (a 11 r)e lip -^P z dtpdr, 
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where D z is defined by 

□ 



if 



d 2 



(7) 



(8) 



• dt 2 dz 2 ' 

The coefficient at the r.h.s. denominators has the follow- 
ing form 

N nt = y [Jf{ai n r ) - Ji-i(ai n r Q )Ji + i(ai n r )] . (9) 

In the case of a single-mode fiber (with two polariza- 
tion) we could omit the difference between group veloci- 
ties, assuming that the fiber is isotropic and do not have 
bending 0- However in a multi-mode fiber it has to be 
taken into account, because different modes have sub- 
stantially different group velocities. We describe it intro- 
ducing the co-ordinate system which moves with average 
group velocity 



t = (t — (3 z)e, 



(10a) 
(10b) 

(10c) 



Next step is introducing a slowly varying amplitude of 
the wave envelope [l II in the form 



-ikz 



Hi) 



and X have the same unit as electrical field [V/m] 



III. TWO MODES INTERACTION 

First, due to the isotropic material assumption, we 
have fcoi = ^oi an d ^jj = &n- Plugging a solution for 
electromagnetic field 0,0] in to the left hand side of the 
equations ©-(0 yields four equations (two modes and 
each mode have two polarizations): 

e 2 k" 

id s X^ - i6d T X± + -^-d TT X± - 



2a 01 

-uj i-^2|^S| 2 + ^31^-111 J ^01 

'or [0 4 \Xl\ 2 + Q 5 X^Xti* + OeX^X**] X±, 



[Oi|^o ±12 



h\x£\ 2 ]x£ 



(12a) 



3 



la 



Pn [Qi|* n | 2 + Q 2 |X n | 2 + Qal^il 2 ] Xu 

(12b) 



where 



^ — o (^01 ^'11) • 



(13) 



The coupling coefficients have the following form 

dr, (14) 
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4fc 2 
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Ji(aur) (Jo(anr) + J 2 2 (aiir)) 



ro 

/ r2J 4 



(aoir)dr, 



dr, 
(15) 

(16) 



ro 

y r2J 1 4 (anr)dr 



2fc 



11 7 2 



3a 2 ! 



Ji (aur) (Jo(anr) + J2(anr)) z 



dr, 
(17) 



+ 



:! = / r2J 1 (aiir)J (a i? , )c?r 
/" 2fc 2 

+ / r^-J 2 (a 01 )(j 2 ( ail r) + J 2 (a n r))dr 
11 

J r^^-J 1 (anir)Ji(a 11 r)Jo(aoir)J 2 (a 11 r)dr 



-J r ^ 01 ^ 11 Ji(ni,i r)./.i('» Ll /')./„ (rim / )./o(<ijj '')<//'. 



3a ian 



(18) 



r2J x (aur) J (a i»")dr 





ro 



4fc 2 

r — -o- J\ («i 1 r) J 2 (aoi 
3a oi 



/ r ^ fc ° lfcn j t (aur) Ji (a ir) J (a ir) J 2 (anr)dr, 
Jaoian 



/ r 2fc ° lfcn Ji(aiir)Ji(aoir)Jo(aoir)J (anr)dr, 



3a iaii 



(19) 



ro 

4 = J r2J?{a n r)J$(a 01 r)dr+ (20) 


/2fc 
r 3af" J ° 2 ^ 01 ^ ( J 2 ( anr ) + J o ( aiir )) dr 



- I r-^™— — J 1 (a ir)Ji(anr)Jo(aoir)Jo(aiir)dr 
oaoian 



^ 2fcoifcn ^ ( aoir ) ^ (aur) Jo(aoir) J 2 (anr)dr, 
3a ian 



ro 

= y r2J 1 2 (a 1 ir)J 2 (a ir) (21) 



4^,2 

Mm 



+ /r- Ji(aiir)Ji(a ir)J (a:oir)Jo(Q!iir)dr 
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ro ^ 

+ I r 01 11 Ji{ai l r)J l {a m r)J {a m r)J 2 {a l ir)dr 1 
3aoiaii 



ro 

h,6 = ©5,6 = y r2J 1 2 (a 11 r)J 2 (a ir) ( ir, (22) 




and 



Pa = 



3w 2 x 2 



32A^ic 2 ' 



(23) 



where I = 0, 1 and P ;i have units [(F 2 m) while x is 
in [m 2 /^ 2 ]. 
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The coefficients Ox and O2 are the same as coefficients 
for one mode fiber with I = and n = 1. The Q x and 
Q2 coefficients for the multi-mode fiber are the same as 
for one mode one with I — 1 and n = 0. This part of 
equations describe interaction between polarization mode 
with same number I and have been analyzed in Q. 

Rest of the coefficients describe the coupling between 
two different modes (with different I number and theirs 
polarizations). 

Notice that coefficient O5, ©6, Q5 and Q>6 have same 
values because they describe mixed interaction (between 
different modes and different polarization). 
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IV. NUMERICAL RESULTS 



FIG. 2: Numerical results for coupling coefficient, difference 
between Q 5 , Q6 and other coefficient. 



First numerically evaluate eigenvalues aoi and an 
from Hondros-Debye equation. In the next steep in- 
tegrals with Bessel function are calculated numerically. 
Eigenvalues are evaluated without approximation of 
weakly guided fiber where we have £\ ~ £2 0- 

The results for inter-mode influence arc shown at the 
figure n Note that the mode 01 have cutoff frequency 
near V ~ 2.4 for that reason the coefficients have dif- 
ferent behaviour. Coefficients O5, ©6, Q5 and Qe de- 



ters with following values: 

lu = 12.2 *10 14 Hz (A « 1.54/im), 

E\ = 2.25 (ref. index 1.5), 

e 2 = 1.96 (ref. index 1.4), 

r from 1.2 * l(T 6 m to 10 * 10~ 6 m. 



CONCLUSION 
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FIG. 1: Numerical results for coupling coefficient for mode 01 
and 11 



scribe mixed interactions (different modes and differ- 
ent polarizations) and for isotropic medium have same 
value. Figure El shows difference between coefficients 
(O5 = 6 = Q5 = Qq is smallest and have smallest incre- 
ment). All of this coefficients are smallest than ©1 — O4 
coefficients. 

In our calculation we defined normalized frequency as 



V = — rWei - £2, 
c 



(24) 



In this paper we considered the influence of Kerr non- 
linearity on the mode coupling in the case of two-mode 
fiber. The results show how the modes influence each 
other and we expect higher influence when V is big- 
ger. The reason is that the propagation constants fcoi 
for higher V is roughly fc n (fig. |3J). 



Core refractive index 




2 4 6 



10 12 14 16 18 20 
V 



and in numerical calculations we used physical parame- 



FIG. 3: Difference in propagation constants between the 11 
mode and the 01 mode, where k — ui/c 

We can reconsider equation l|12|l again but with sim- 
plification, that we do not take into the account polar- 
ization interaction. Choosing two modes with the same 
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polarization yield: 

e 2 k" 

id^X Q i - i5d T X Q i + — — d TT X Q i = 

ZO~ 

Poi [Oi|X i| 2 + O 3 |Xii| 2 ] Xoi, (26a) 
e 2 k" 

ideXn + i8d T Xu + — — d TT Xu = 

2(7 

Pii [Qil^nl 2 + Qsl^oil 2 ] Xu. (26b) 

In this case we have only four parameters (two for 
self phase modulation and two for cross phase modula- 
tion), which describe only interaction between modes. In 
boundary when we have two modes with the same group 
velocity (5 = 0), we get system of equations like this in 
case of polarization interaction |6|. 



Method which we used for a cylindrical waveguide is 
easily reformulated for a different waveguide shape where 
second-modes (and higher modes) are important, for ex- 
ample in the elliptical waveguide p|. This method also 
can be used for a photonic crystal fiber , dir ectly or within 
the approximations of the approach of [Tjj . Here we have 
delivered all calculations for isotropic medium but there 
is a possibility to make it for anisotropic medium. 
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